Critical Behavior and Extended States in 2D and 3D Systems with Gas-like Disorder 
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With a tight binding treatment we examine amorphous conductors with gas-like disorder, or no 
correlations among the site positions. We consider an exponentially decaying hopping integral with 
range I, and the Inverse Participation Ratio (IPR) is used to characterize carrier wave functions 
with respect to localization. With the aid of two complementary finite size scaling techniques 
to extrapolate to the bulk limit (both methods exploit critical behavior in different ways to find 
the boundary between domains of extended and localized wave functions) which nevertheless yield 
identical results, we obtain phase diagrams showing regions where states are extended and domains 
of localized states. In the 2D case, states are localized below a threshold length scale l c on the order 
of the interparticle separation p _1,/2 with a finite fraction of states extended for / > l c . For D = 3, 
the extended phase is flanked by regions of localized states and bounded by two mobility edges. The 
swath of extended states, broad for I ~ 1, becomes narrower with decreasing I, though persisting 
with finite width even for I < ip -1 / 3 . Mobility edges are interpreted as lines of critical points, and 
we calculate the corresponding critical exponents. 
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PACS numbers: 72.15.Rn, 72.80.Ng, 71.23.-k, 71.23.An 

The Anderson model with disorder manifested as 
random perturbations in the site potential (i.e. diago- 
nal disorder) has provided a way to examine localiza- 
tion effects in the context of a tight binding treatment. 
However, off-diagonal disorder where there are stochas- 
tic variations in the tunneling matrix elements also is 
germane to situations examined in experiment [2[ and 
theoretical studies 0-0 which have also included calcu- 
lations of the energy density of states (DOS) Iq-|l0l| and 
aspects of localization for 3D geometries [lll - ll5j . Our 
aim is to determine conditions for the existence of ex- 
tended states (system spanning wave functions amenable 
to charge transport) in the presence of very strong dis- 
order where there is no correlation among site positions, 
a circumstance known as gas-like disorder [l6| and also 
termed topological disorder due to the absence of a spe- 
cific underlying lattice geometry. 

Setting the constant diagonal terms to zero, we con- 
sider the tight binding Hamiltonian 
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ranges over the N 
particles contained in the simulation volume, we take 
the hopping parameter to to be 3.0 electron volts, and 
the factor of "1/2" compensates for multiple counting. 
The operators c' and c create and destroy occupied elec- 
tronic orbitals at sites indicated by the subscript. For 
the hopping integral, we use an exponential dependence 



where ru is the separation between 



sites i and j, s = p~ x l D is the typical distance between 
neighboring sites, p is the volume density of sites, D is 
the dimensionality of the system, and 7 is a dimension- 
less parameter. Since the length scale for the decay of the 
tunneling matrix element is I = s/7, large/small 7 val- 



ues correspond to decay lengths small/large in relation to 
the typical interatomic separation. For the sake of con- 
venience, we rescale coordinates such that p = 1, with a 
simple inverse relationship I — 7 , and V(ry) — e~ JTij . 
Although the matrix element V(ry) is finite in range by 
virtue of the exponential decay, we nevertheless take into 
consideration hopping among all pairs of orbitals con- 
tained in the system of N sites with only a negligible 
increase in computational burden (i.e. a contribution on 
the order of A^ 2 ) relative to direct diagonalization used 
to calculate carrier states, which scales as A^ 3 . 

With positional order (either locally or globally) not 
preserved in gas-like amorphous materials, the disorder 
strength is nonetheless characterized in a sense by the 
hopping integral decay scale I. For large I (small 7), 
the comparatively long-ranged tunneling connects sites 
to many neighbors, effectively averaging over hopping 
rates to and from many sites and mitigating the effect 
of strong disorder. However, if I <g; 1 (large 7), tunneling 
chiefly involves nearest neighbors. Hence, disorder is an 
a sense amplified for 7 > 1 and muted if 7 <C 1, and 
one proceeds with the intuition that extended states are 
more likely to exist for larger I than for smaller I values. 

For 2D systems, we find a threshold l c where all states 
are localized for I > l c ~ 1, with a finite portion of the 
wave functions extended for I < l c . In the 3D case, a 
meta-study Q of a large body of experiments in dis- 
parate systems corresponding to a broad range of 7 val- 
ues suggests a termination of conducting behavior for 
7 > 3.8. Nevertheless, though we find the portion of 
extended states to diminish monotonically in 7, we do 
not find an abrupt termination of the extended region 
for any particular value of 7, in contrast to D — 2 where 
extended states do not exist if I is smaller than the inter- 
particle separation. We construct phase diagrams show- 



2 



ing a region of extended states for intermediate energies 
bounded by two mobility edges; for a specific 7 value, the 
latter are critical points marking a second order transi- 
tion from the localized to extended phase at the lower 
edge and from extended to localized states at the upper 
edge. 

To implement the numerical calculations, we examine a 
L D supercell using periodic boundary conditions to miti- 
gate finite size effects. A range of L values (i.e. as high as 
L = 110 for D = 2 and L = 22 for D = 3) are used to per- 
form finite size scaling to determine if electronic states are 
localized or extended in the bulk limit. The Inverse Par- 
ticipation Ratio, Y 2 = J2iLi lV>i| 4 /(2iLi l^il 2 ) 2 shows 
distinct behavior depending on whether the wave func- 
tion is confined to a small volume (larger Y 2 ) or spread 
out over a larger region (smaller Y 2 ), and hence more ex- 
tended in character. With a crystal lattice absent, one 
operates in the grand canonical ensemble, and the num- 
ber of sites TV varies about (N) = pL D from sample to 
sample. In generating realizations of disorder, we account 
for fluctuations in the numbers of particles in an unbiased 
way using a technique described elsewhere 17 1. To ob- 
tain adequate statistics, 10 5 wave functions are retained 
for each combination of L and 7. 

The random character of disorder precludes the study 
of the evolution of individual states with increasing L, 
and instead aggregates of wave functions must be exam- 
ined. With carrier states parameterized by energy eigen- 
values, one could in principle create partitions of width 
8E centered about uniformly spaced energies. However, 
since the energy density of states is sharply peaked fl7^ . 
energy channels far from the central peak suffer from 
poor statistics. We ensure uniform statistics by using 
normalized rank f = r/M (ranging from to 1) to la- 
bel carrier states, where M is the total number of states 
and r is the global eigenvalue rank within the large aggre- 
gate. For n = 100 partitions of the f domain, the channel 
width Sf is sufficiently small that systematic admixture 
effects from neighboring channels play a minor role, while 
still providing sufficient statistics for analysis, and paral- 
lel calculations for n = 50 and n = 200 yield results in 
quantitative agreement with the n = 100 scheme. 

We calculate the channel averaged IPR in the ther- 
modynamic limit as an avenue for the identification of 
mobility edges. For moderate to large L, we use a power 
law formula Y 2 {L) = Y§ + A X L'^ + A 2 L- & where K>° 
is the participation ratio in the bulk limit, /3 and S are 
leading and next-to-leading order exponents, and Ai and 
A 2 are amplitudes. The parameters Y 2 °, /3, 5, Ai, and A 2 
are fixed by a Monte Carlo nonlinear least squares fit to 
the channel averaged participation ratios. 

Extrapolated participation ratios appear in Fig. [1] The 
troughs in the bulk Y 2 curves shown in panel (a) of Fig. Q] 
deepen with decreasing 7; while IPR minima for 7 > 1.0 
are manifestly nonzero, participation ratios do seem to 
descend to zero for less rapid decay rates. Neverthe- 
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FIG. 1: (Color Online) Extrapolated Participation Ratio 2D 
and 3D results appear in panel (a) and panel (b), respectively 
for various 7 values. 



less, the gradual approach to the abscissa confounds a 
straightforward objective identification of a mobility edge 
for the 2D systems we examine. 

Irrespective of 7, the D = 3 extrapolated IPR curves 
in panel (b) Fig.Q]seem to plunge to zero and remain zero 
for a range of f values. However, the width of the interval 
where extended states seem to exist decreases with 7, 
with Y 2 appearing to vanish for most of the f domain in 
the case 7 = 1.0 while only briefly touching the abscissa 
for 7 = 5.5. For 7 < 3 the IPR falls to zero with no 
subsequent return to finite values. On the other hand, 
for 7 > 3, the interval of vanishing participation ratios is 
terminated for r < 1 as Y 2 abruptly rises to finite values. 
Salient feature not present in the 2D results are derivative 
discontinuities where Y 2 falls to zero, descending to or 
rising from the horizontal axis with a finite slope. The 
clear delineation of the regions where the IPR vanishes 
favors the use of basins in the Y 2 curves as a way to 
identify the boundary between extended and localized 
states in the 3D case. 

With a complementary analysis, we calculate phase 
portraits showing localized and extended states for D = 2 
and D = 3. As a consequence of critical behavior at the 
mobility edge, one expects from single parameter finite 
size scaling theory a form Y 2 — L a ° FulL 1 ^ (f — r c )] 
near the boundary between regions of extended and 
localized states, with an and vd being critical expo- 
nents (vd controls the singular behavior of the corre- 
lation length £). Instead of trying to determine numer- 
ically where Y 2 falls to zero, the mobility edge is iden- 
tified in an objective fashion by locating intersections of 
<f) 2 = Y 2 (f, L) _1 L _QD , where an is fixed by insisting that 
curves for different system sizes L coincide, with f c deter- 
mined by the location of the crossing. A similar approach 
has been applied in theoretical studies of a 2D Anderson 
model with long-range correlated disorder 18[ and An- 
derson localization of phonons 19] Data collapses, where 
4> 2 points plotted versus L 1 ^ (f — f c ) for various L coin- 
cide on a single curve also are associated with a critical 
point, and in the case of the lower mobility edge are sharp 
enough to permit the determination the critical index vo 
by optimizing the data collapse. 

Across the range of 7 examined, we find from the 4> 2 
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FIG. 2: (Color Online) Data collapse plots for (top to bottom) 
7 = 0.3, 7 = 0.5, and 7 = 0.75 appear on the left with 
the corresponding (j>2 intersections displayed on the right for 
D = 2. 



curve crossings that a 2 d — 1.6 ± 0.1 for D = 2 and 
c*3d = 1.25 ± 0.15 for D = 3. In the 2D case, a mo- 
bility edge is identified only for 7 < 1, and in the right 
panels of Fig. [21 </>2 intersections for four distinct L val- 
ues are shown for 7 = {0.3, 0.5, 0.75}. The left panels 
show good data collapses for ^2d(7 = 0.3) = 2.3 ± 0.2, 
^2d(7 = 0.5) = 2.2 ± 0.2, and */ 2 d(7 = 0.75) = 2.3 ± 0.2 
The i>2d exponents in each case are identical within the 
bounds of error, in accord with the Harris Criterion. Al- 
though only one intersection appears for 7 = 0.3 and 
7 = 0.5, the (j)2 curves appear to intersect twice for 
7 = 0.75 with the second crossing near f = 1.0. The 
incipient upper mobility edge is mirrored in the data col- 
lapse, which for 7 = 0.75 is defocused in the vicinity 
of f = 1.0, where the scaling of Y2 is controlled by the 
critical point marking the upper mobility edge. 

A similar treatment identifies mobility boundaries for 
D = 3, and the corresponding 4>2 data collapses and in- 
tersections are shown in Fig. [3] For 7 < 2.5, only a single 
intersection marks the boundary between localized states 
for r < r c and extended states for f > f c . On the other 
hand, two sets of (f>2 intersections may be discerned for 
7 > 2.5, indicating two distinct mobility edges with an 
interval of extended states flanked by localized states for 
f < f l ° wor and f > f" ppcr . In the context of a similar- 
tight bindin g m odel, a calculation by J. J. Krich and 
A. A-Guzik [14| has also identified two sets of mobility 
edges for D = 3. The case 7 = 2.5 is marginal, since the 
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FIG. 3: (Color Online) Data collapse plots for various 7 values 
(7 increases from top to bottom); data collapse plots corre- 
sponding to the lower mobility edge appear on the left, with 
4>2 intersection plots shown on the right for D — 3. 



upper intersection appears to coincide with the upper f 
extreme, r — 1. The convergence of the mobility edges, 
and the concomitant constriction of the extended region 
affects the upper part of the data collapses calculated for 
the lower mobility edges where the severity of the defo- 
cusing increases with the proximity of the critical points 
in the upper mobility boundary. The areas in the right 
panels marked in gray match the f domains between in- 
tersections of <j>2 curves, and are determined from the 
extrapolated IPR results in Fig. [TJ where we adopt the 
criterion Y2 < 10~ 3 for the presence of extended states. 
For the lower mobility edge, scaling collapses are sharp 
enough to be of service in fixing the exponent i/3d, with 
results shown in Table |H As in the 2D case, within error 
bounds, the critical indices are in agreement across the 
broad range of 7 values under consideration. 

Phase diagrams showing regions of extended and local- 
ized states for D = 2 appear in Fig. 21 where the main 
graph shows results in terms of the normalized rank f, 
and the inset is a phase portrait expressed in terms of the 
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TABLE I: results with for various 7 values 
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FIG. 4: Phase portraits for D = 2 in terms of the normalized 
channel number f. The inset phase diagram is calculated with 
respect to energy. Filled symbols are calculated, and broken 
lines are extrapolations. 



carrier state energies. Extended states occupy the upper 
left corner of the phase diagram, for 7 < 1.0 and in the 
upper E and f ranges. 

Phase portraits rendered in terms of f and E for the 
3D case appear in the upper half of Fig. [5] mobility edges 
calculated from §2 intersections and gleaned from bulk 
Y% values appear in black and gray, respectively; the close 
overlap is a sign of excellent agreement among results of 
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FIG. 5: Phase portraits for D = 3 are shown in the upper 
row; mobility edges obtained from scaled IPR curve inter- 
sections and extrapolated participation ratios appear on the 
same graph. The logarithm of the width of the extended zone 
is graphed with respect to 7 in the lower panels 



the two very distinct methods of locating mobility edges, 
and serves to validate the use of bulk Y 2 results in con- 
structing phase diagrams for 3D systems. The lower part 
of Fig. [5] contains graphs of the logarithm of the widths 
(7) and we (7) of the interval of extended states. Al- 
though the extended region grows rapidly narrower for 
7 > 3.0, there is no indication of an abrupt termination 
for finite I; the widths, within the bounds of error, both 
are consistent with an asymptotically exponential scaling 

e -A/l f or I ^ 1 

In conclusion, we have obtained phase diagrams for 
2D and 3D amorphous systems with gas-like disorder, 
finding extended states to be supported in the upper part 
of the energy range for D — 2 above a threshold range 
l c ~ 0.9 with the appearance of a second mobility edge 
near l c . On the other hand, for D = 3, we find a region 
of extended states even for I <C 1 with no evidence for a 
termination of the extended phase for finite I. Mobility 
edges are interpreted as lines of critical points, which 
we have characterized by calculating the corresponding 
critical exponents. 
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